Correlation between a-decay Energies of Superheavy Nuclei Involving Effect of 

Symmetry Energy 
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A formula for the relationship between the a-decay energies (Q values) of superheavy nuclei (SHN) 
is presented, which is composed of the effects of Coulomb energy and symmetry energy. It can be 
employed not only to validate the experimental observations and measurements to a large extent, 
but also to predict the Q values of heaviest SHN with a high accuracy generally which will be very 
useful for future experiments. Furthermore, the shell closures in superheavy region and the effect 
of the symmetry energy on the stability of SIfN against a-decay are discussed with the help of this 
formula. 
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The synthesis and identification of superheavy nu- 
clei (SHN) have been receiving a worldwide attention 
since the prediction of the existence of superheavy is- 
land in 1960s. But where the closed shells are located 
in the superheavy region is less certain, depending on 
the model employed. The experimental investigations 
are thus crucial and a series of experimental efforts so 
far have been focused on the direct production of SHN 
in heavy ion fusion reactions. The superheavy elements 
with Z ~ 107— 112 have been successfully produced at 
GSI, Darmstadt, in cold-fusion reactions [ij. Several new 
elements with Z = 113—118 have been discovered at 
JINR-FLNR, Dubna, using hot-fusion evaporation reac- 
tions with the neutron-rich ^®Ca beam and actinide tar- 
gets The element 114 was independently confirmed 
recently by the LBNL in the USA [| and GSI Q. A su- 
perheavy element isotope ^^^114 was observed in LBNL 
last year [s], and an isot^e of Z = 113 has been iden- 
tified at RIKEN, Japan Q. Thus up to now superheavy 
elements with Z = 104—118 have been synthesized in ex- 
periment and hence it offers the possibility to study the 
heaviest known nuclear island of stability with greater 
detail. 

Superheavy nuclei with atomic numbers beyond 110 
predominantly undergo sequential a-decay terminated by 
spontaneous fission T| , leading to a-decay that is one ef- 
ficient approach to identify new nucleus via the obser- 
vation of a-decay chain and to extract some information 
about their stability. In experiment one usually measures 
the a-decay Q values and half-lives, while one of the ma- 
jor goals of theory is to predict the half-lives to serve 
the experimental design. As one of the crucial quantity 
for a quantitative prediction of decay half-life, Q value 
strongly affects the calculation of the half-life due to the 
exponential law. Therefore, it is extremely important 
and necessary to obtain an accurate theoretical Q value 



in a reliable half-life prediction during the experiment 
design. However, the existing microscopic nuclear many- 
body approaches do not achieve a very good accuracy. In 
this study, we propose a new approach to calculate the 
a-decay energy with a high accuracy for the superheavy 
elements above 110. In our previous work Q, a formula 
was proposed for a-decay Q value of SHN based on a 
liquid drop model. Taking no account of the shell energy 
it gives as 
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with a = 4ac/3 = 0.9373 
e = -27.4530 [7]. Here 



-4a., 



(1) 



-99.3027 and 



b 

Z, N and A are the proton, 
neutron and mass numbers of the parent nuclei, respec- 
tively. The first two terms on the right hand side are the 
contributions of Coulomb energy and symmetry energy, 
respectively. The nuclear symmetry energy plays an im- 
portant role in astrophysics [8|, |9j] , the structure of exotic 
nuclei and the dynamics of heavy ion reactions (lol - [l^ . 
In this Letter, the effect of symmetry energy on the sta- 
bility of SHN against a-decay is going to be shown. 

Here we study the relationship between the Q values of 
the neighboring SHN taking Eq. ^ as the starting point 
but we do not use the parameters in Ref. any longer. 
With /? = {N — Z)/A denoting the isospin asymmetry 
and Z = A{1 - fi)/2, we obtain 
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Once the decay energy Qi of a reference nucleus "^Zi is 
known, the Q2 values of the other nucleus ^Z2 (target 
nucleus) with the same mass number A can be estimated 
by 



{P2-P1) 



-a,A^'^(P + 2) + ^a,y^P 
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with P — [Pi + /?2)/2 and ac=0.71. The mass depen- 
dence of the symmetry energy coefficient is given by 
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Danielewicz and Lee |l3j | as Osym — Cgym 
where Cgym is the volume symmetry energy coefficient 
of the nuclei and k is the ratio of the surface symme- 
try coefficient to the volume symmetry coefficient. Here 
Cgym = 31.1 and k = 2.31 are taken from the results of 
Rcf. [14] without including the uncertainty. 

Apart from discussed above, (^,/3) or {N,(3) can 

be also adopted as variables. By an analogous deriva- 
tion, the correlation between the Q values of the nuclei 
belonging to an isotope chain with a proton number Z is 
given by 



^a,Z^/''{l - /3)-2/3(i + 2/3) + 8a,y„,/3 



(4) 



and that of the nuclei belonging to an isotone chain with 
a neutron number TV is given by 

Q2 - Oi - {132 ~ Pi) X 
25/3 

— acN^'^l + /3)-5/3(ll + 5^ + 2/32) + Sasyra^ 

y 

In general, if one selects ^ = xZ + yN and /3 as variables, 
the relationship between the Q values of a-decay can be 
written as 



Q2 = Qi 



- {P2 - /3i 



X {^«cf /^[(l 



P)2 
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(ll + 5/3 + 2/32)y] +8a,j,™/3}. 
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where x and y are integers and |x| + 7^ with 
Z = {I - (i)^/[{l-(3)x + {l+(3)y] and iV = (1 + 
[(1 - P)x + (1 + P)y]. Here only the differences of 
the symmetry energy effect {asym term) together with 
the differences of Coulomb energy effect {uc term) be- 
tween a reference nucleus and a target one contribute to 
this correlation. The isospin dependence of the symme- 
try energy coefficient agym is neglected here because the 
asym changes quite slightly between the neighboring nu- 
clei. With this formula, the Q values of target nuclei can 
be obtained by any neighboring nuclei. Eqs. (01) and 
^ are the special cases of Eq. (|6])-that is, x — y = 1 for 
Eq. (P, a; = 1, 2/ = for Eq. ^ and x = 0, y = 1 for 
Eq. (O, respectively. 

In order to test the applicability of Eq. ([5]), we com- 
pute the Q values of recently synthesized heaviest SHN 
with the help of their neighbors, and the results are listed 
in Fig. 1 compared with experimental ones. The results 
obtained with Eqs. ([3]), (|1]) and ([S]) are marked by distin- 
guishable symbols. For the nuclei except '^^'^118, ^^"llS, 

113 and ^^^111, our approach reproduces the mea- 
sured values quite accurately with a root-mean-square 
deviation ^/(a^ = 0.077 MeV and an average deviation 
(a) = 0.064 MeV for central values from 380 reference- 
target combinations. It is thus very practical that Eq. 
(|n]) can be reliably applied to the Q values of the as-yet- 
unobserved SHN with the help of known nuclei which is 



the most effective method to the Q values at present. As 
three simple cases of Eq. ([B]), Eqs. ^ and ([5|) work 
even better with y/Ja^ = 0.052 MeV and (a) ^ 0.043 
MeV from 96 reference-target combinations, which are 
very convenient to be used and are sufficient for predic- 
tions of Q values generally, though they are simple in 
formalism. In addition, the agreement between the ex- 
perimental and theoretical values has additional signif- 
icance. Since the Q values of the reference nuclei are 
taken from the experimental measurements in calcula- 
tions, the agreement suggests that the experimental data 
themselves are consistent with each other, which indi- 
cates that the experimental observations and measure- 
ments of the SHN are reliable to a great extent. These 
SHN still await independent verification by other labo- 
ratories, which is not easy because the new SHN form 
an isolated island that tends to be not linked through 
a-decay chains with any known nuclei, making the the- 
oretical supports become important and necessary. For 
the nuclide ^^"llS, the experimental value of 10.14±0.41 
j-j MeV carries a large uncertainty while the Q value is 
about 10.4 MeV according to Eq. ([5]), which requires 
a more precise experimental measurement. 

The shell closures should play a particular important 
role in the superheavy system. However, modern theo- 
retical approaches disagree on the position of the closed 
shells. For instance, the macroscopic-microscopic models 
with various parameterizatio ns p redict the shell gaps at 
Z = 114 and N = 184 [H [ll. Skyrme-Hartree-Fock 
calculations favor Z = 124, 126 and N = 184 [13, d 
while the relativistic mean field models favor Z = 120, 
N = 172 [il-i^ and Z = 120, N = 184^. The magic 
numbers Z = 132 and N = 194 were predicted from the 
discontinuity of the volume integral at shell closures [2^ . 
The reason for this uncertainty lies in incomplete knowl- 
edge of the nuclear force and the difficulty of many-body 
techniques. It is well known that the shell effect on the a 
radioactivity is related to the Q value. For the a-decay 
of the nuclei being not close to the shell closures, due 
to a parent nucleus and its daughter nucleus sharing the 
same odevity of both the proton and neutron numbers, 
the shell correction (also pairing correction) energies to 
their masses could be canceled to a large extent leading 
to a small correction to a Q value compared with the con- 
tributions of the Coulomb and symmetry energies within 
semi-empirical formulas Q, and even these small shell 
energies to the Q values turn out to be nearly a constant 
in a local region of particle numbers which confirms that 
the shell energies hardly take effect in Eq. ([5]). Most 
importantly, the agreements between the estimated and 
experimental results in turn show this point. Once the 
parent nucleus or the daughter one has neutron and/or 
proton magic numbers or the shell gaps are crossed, the 
Q value shows an irregular behaviour. Since the shell 
energy is excluded in Eq. (|6]), it should show some dis- 
crepancies for nuclei around shell closures. Yet, it would 
help us to investigate the shell structure by comparing 
the experimental and calculated Q values. All the theo- 
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FIG. 1: Comparison of a-decay Q values with Eq. ((6]) (the rectangles with error bars) and experimental ones 01 (shaded area) 
of recently synthesized heaviest SHN. The horizontal ordinate denotes the mass numbers of the reference nuclei. Since the 
experimental Q values of the reference nuclei include uncertainties, the calculated ones also display error bars. The results from 
Eqs. ([3}, (|4]) and ([5)l which are special cases of Eq. ([6}, are presented separately marked by hollow rectangles. 



retical calculated Q values of '^^^11^, 282113 and ^som 
based on Eq. ^ are lower than the experimental ones, 
which is possibly attributed to the likely shell gaps at 
Z = 120 for 294118^ and at = 166 for 282113 and 
28"111. In Ref. [23], it is suggested that N = 166 is a neu- 
tron shell gap in a certain region within relativistic mean 
field models. Yet, to confirm the existence of shell gaps 
positively is not easy due to the insufficient experimen- 
tal observations. Apart from the shell effects, dramatic 
shape changes could also affect the energy of a-decay, as 
pointed out in Ref. [23|. Nevertheless, the deviations are 
not more than 0.5 MeV in general. However, it should 
be much easier to confirm the non-existence of shell clo- 
sures. Once the Q values behave quite regularly in a 
local range, a magic number should not appear here. For 



the eight nuclides of elements 116 and 114 (2^" 293^15 
and 286-289 ^^-j^^-j together with the six nuclei with a neu- 
tron number N = 174 (290ll6, 289115 and 288114) and 
N = 172 (28^115, 286^14 a^j^d '^s.^Wi)^ the experimental 
Q values can be reproduced very accurately that con- 
firms Z — 114 and = 172 are not shell closures in the 
considered region. Of course, one cannot rule out the 
possibility that they appear as magic numbers in other 
mass regions. 

We now turn to the effect of the symmetry energy on 
a-decay. The a^ym term (difference of the symmetry en- 
ergy effect between the reference and target nuclei) con- 
tributes by about 45% in Eq. (P, 80% in Eq. (jj) and 
35% in Eq. (O to the AQ = Q2 - Qi for the SHN in 
Fig. 1, suggesting its important role in the correlations 



4 



-1 — ' — I — ' — I — ' — I — ' — r 
0.0k^„ Z=114 




J 1 I 1 L 



-1 — ' — r 
0.0 <v 



-I 1 1 1 r- 

Z=116 



— *— Exp.\ <»x^ 
— o — « term\^ <* 

— o — a term \ 

sym 

— • — a + a term 



286 288 290 292 294 
mass number 



290 292 294 296 298 
mass number 



Qi in Eq. Q taking the elements 114 and 116 as examples 
(^*®114 and ^^'^116 as reference nuclei with decay energies Qi, 
respectively). The experimental data, if available, are also 
shown for comparison. 



between the Q values of SHN. The large contribution 
of Ugym term in Eq. @ is of particular importance as 
discussed below. The most significant experimental con- 
clusion is these observed superheavy elements generally 
display a trend of increased stability with larger neutron 
number, which is almost attributed to the larger symme- 
try energy that lowers the Q values. In order to illumi- 
nate this conclusion more obviously, we plot in Fig. 2 
the contributions of the a^ym and Uc terms to Q2 — Qi 
in Eq. One can find that the agym term contributes 
much more greatly than the Uc term. Therefore, due 
to the inclusion of the effect of symmetry energy, apart 
from the theoretical estimations being able to agree with 
the experimental values, the Q values reduce much more 
rapidly as N increases, and hence a superheavy element 



becomes longer-lived against a-decay with increasing N. 
In other words, it is the symmetry energy that primarily 
enhances the stability against a-decay with larger neu- 
tron number for these synthesized SHN not around shell 
closures. 

We have investigated some aspects of the a-decay of 
SHN. The main conclusions are summarized as follows: 
(1) A simple formula for the correlation between the a- 
decay Q values of the SHN has been proposed, which 
works very well for an estimation of the a-decay energies 
of the recently synthesized SHN. They thus allow us to 
reliably predict the Q values of the still unknown SHN 
with a good accuracy, and is going to be very useful for 
the future experiment design. Also, the agreements be- 
tween the calculated and experimental values indicate the 
reliability of the experimental observations and measure- 
ments on these synthesized SHN to a great extent. (2) 
Z = 114 and N = 172 turn out to be not shell closures 
for the presently observed superheavy region experimen- 
tally. (3) The observed increase of a-decay half-lives with 
increasing neutron number, i.e., the increased stability of 
these SHN not around shell closures with larger neutron 
number, is primarily attributed to the effect of the sym- 
metry energy. 
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